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Abstract 

Frustration of magnetic systems which is caused by competing interactions is the 
driving force of several unusual phenomena such as plateaus and jumps of the mag- 
netization curve as well as of unusual energy spectra with for instance many singlet 
levels below the first triplet state. The antiferromagnetic cuboctahedron can serve 
as a paradigmatic example of certain frustrated antiferromagnets. In addition it has 
the advantage that its complete energy spectrum can be obtained up to individual 
spin quantum numbers of s = 3/2 (16,777,216 states). 
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1 Introduction 



The magnetism of antiferromagnetically coupled and geometrically frustrated 
spin systems is a fascinating subject due to the richness of phenomena that are 
observed [T1I2] . Realizations of such systems exist in one, two, and three dimen- 
sions; the most famous being the two-dimensional kagome lattice [2,3,4'' ^|5f7] 
and the three-dimensional pyrochlore antiferromagnet pPllUfll|[T2iil3^14il5lll6] . 

It is very interesting and from the point of theoretical modeling appealing 
that similar but zero-dimensional spin systems - in the form of magnetic 
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Fig. 1. Cuboctahedron, icosidodecahedron, and (part of the) kagome lattice 




molecules [T7|18II19II20|21] - exist that potentially could show many of the 

special features of geometrically frustrated antiferromagnets. Figure [1] dis- 
plays the zero-dimensional "little brothers" of the kagome antiferromagnet: 
the cuboctahedron which consists of squares surrounded by triangles and the 
icosidodecahedron which consists of pentagons surrounded by triangles. Such 
finite size antiferromagnets offer the possibility to discover and understand 
properties that are shared by the infinitely extended lattices. An example is 
the discovery of localized independent magnons [61122] , which explain the un- 
usual magnetization jump at the saturation field. Also the plateau at 1/3 of 
the saturation magnetization that appears in systems built of corner sharing 
triangles could be more deeply investigated by looking at the cuboctahedron 
and the icosidodecahedron [23)121] . 

In this article we continue investigations along this line. We focus on two 
points. First we discuss the physics of the regular cuboctahedron as a func- 
tion of the single spin quantum number s = 1/2,1,3/2. For these cases all 
energy eigenvalues could be obtained with the help of Irreducible Tensor Op- 
erator (ITO) techniques [2^26f27j and by application of point group sym- 
metries. As a second point we investigate irregular cuboctahedra. This study 
is motivated by recent magnetization measurements of the icosidodecahedral 
molecules {Mo72Fe3o} [H] and {Mo72Cr3o} [20] published in Ref. [28] which 
could successfully be interpreted by a classical Heisenberg model with random 
antiferromagnetic exchange couplings between the paramagnetic ions. 



2 Theoretical model 

The physics of many of the mentioned spin systems can be well understood 
with the help of the isotropic Heisenberg model, 
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Here the sum runs over pairs of spins given by spin operators s at sites u and 
V. A negative value of the exchange interaction J„„ corresponds to antiferro- 
magnetic couphng. We refer to a regular body, e.g. cuboctahedron, if there are 
only nearest neighbor couplings of constant size J. In the case of an irregular 
couphng the nearest-neighbor couplings can assume values according to the 
chosen distribution. 

Since the Hamiltonian commutes with the total spin, we can find a common 
eigenbasis { | z/)} of 5^, and and denote the related eigenvalues by E^, 
S^, and Mj,, respectively. The eigenvalues of ([T]) are evaluated in mutually 
orthogonal subspaces TCi^S, M) of total spin S and total magnetic quantum 
number M using Irreducible Tensor Operator (ITO) techniques [25|26ll27j . In 
addition point group symmetries have been applied for the regular cubocta- 
hedron. 



3 The regular cuboctahedron 

The regular cuboctahedron belongs to the class of geometrically frustrated 
antiferromagnets built of corner- sharing triangles. Such systems possess an 
extended magnetization plateau at 1/3 of the saturation magnetization A^sat 
caused by dominant up-up-down contributions [23112^ . an unusually high jump 
of the magnetization at the saturation field due to independent magnons [22|6] 
as well as low-lying singlets below the first triplet level [29|30II24] . These fea- 
tures are shared for instance by the icosidodecahedron and by the kagome 
lattice. 
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Fig. 2. Magnetization as a function of applied field at T = for the regular cuboc- 
tahedron with s = l/2, s = l, and s = 3/2. The extended plateau at A4sat/3 is 
clearly visible. 

Figure [2] shows the magnetization curves at T = for the regular cubocta- 
hedron with s = 1/2, s = 1, and s = 3/2. These curves show besides the 
plateau at A^sat/3 a jump to saturation of height AM = 2. Both features 
are reflected by the differential susceptibility function which is displayed in 
Fig. [31 Each step in Fig. [2] corresponds to a peak in Fig. [31 One notices that 
the peaks are washed out for higher temperatures, but that the minimum that 
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corresponds to the plateau at A^sat/3 persists up to temperatures of the order 
of the exchange couphng. 
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Fig. 3. Differential susceptibility as a function of applied field at 
ksT l\J\ = 0.25,0.5,0.75,1.0 for the regular cuboctahedron with s = 1/2, 
s = 1, and s = 3/2. The smoother the curve, the higher the temperature. 

As a function of the intrinsic spin s the differential susceptibility dA^/d-B ex- 
hibits two properties. With increasing spin quantum number s the individual 
peaks oscillate more and more with smaller relative amplitude, but the mini- 
mum at 1/3 is actually sharpened. It is known that in the classical limit, i.e. for 
s — s> oo, the differential susceptibility is practically fiat below the saturation 
field except for the dip at 1/3 [23J. 
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Fig. 4. Low-lying energy levels for the regular cuboctahedron with s = 1/2, s = 1, 
and s = 3/2. Numbers attached to selected levels denote their multiplicities ds] 
unlabeled levels below the highest labeled level have ds = 1. 

For zero field Fig. H] shows the low- lying energy levels. In the case of s = 1/2 
(l.h.s. of Fig.Hj) one notices the low-lying singlets below the first triplet. These 
states are a cornerstone of geometric frustration and as well present in the 
kagome lattice and the icosidodecahedron with s = 1/2 [2U]. It is interesting to 
note that with increasing s, i.e. towards a more classical behavior, the number 
of these states decreases. For s = 1 (middle of Fig. Hj) the first excited singlet 
level is already (slightly) above the lowest triplet level. For s = 3/2 (r.h.s. 
of Fig. H]) a doubly degenerate excited singlet level remains below the lowest 
triplet, the others have disappeared. This behavior, i.e. no excited singlets be- 
low the lowest triplet for integer spins and a doubly degenerate excited singlet 
below the lowest triplet, does not change anymore for higher spin quantum 
numbers as can be checked e.g. by Lanczos methods. 

The rather high symmetry of the cuboctahedron leads to many degenerate 
energy levels. As examples we label some low-lying energy levels in Fig. H] 
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by their multiplicity ds, i.e. by the degeneracy of the whole multiplet. The 
full degeneracy including the multiplicity of the magnetic sublevels (Im is then 
d = dsX dM- Clearly, such high multiplicities have an important impact on the 
magnetocaloric behavior since they increase the entropy for low temperatures 
[30f3T] . In the following we would like to discuss the impact of low- lying singlets 
below the first triplet which in the case of extended lattices are supposed to 
condense in infinite number onto the ground state. 
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Fig. 5. Heat capacity (l.h.s.) and zero-field susceptibility (r.h.s.) for the regular 
cuboctahedron with s = 1/2, s = 1, and s = 3/2. 

Figure [5] compares the heat capacity (l.h.s.) and the zero-field susceptibility 
(r.h.s.) for the regular cuboctahedron with s = l/2,s = l, and s = 3/2. The 
heat capacity shows a pronounced double peak structure for s = 1/2 and s = 1 
which dissolves into a broad peak with increasing spin quantum number. The 
broad peak also moves to higher temperatures with increasing s. The reason 
for the first sharp peak is twofold. Since there are several gaps between the 
low-lying levels the density of states has a very discontinuous structure which 
results in the double peak structure. For s = 1/2 the low- lying singlets provide 
a very low-lying non-magnetic density of states which is responsible for the 
fact that the first sharp peak is at such low temperatures. For s = 1 the first 
sharp peak results from both excited singlet as well as lowest triplet levels. 
For s = 3/2 a remnant of the first sharp peak is still visible; it is given by 
the low-lying singlets, but since they are so few, also influenced by the lowest 
triplet levels. 

The behavior of the heat capacity is contrasted by the susceptibility on the 
r.h.s. of Fig. [5] which reflects mostly the density of states of magnetic levels 
and is only weakly influenced by low-lying singlets. Therefore, the first sharp 
peak, or any other structure at very low temperatures, is absent . 



4 The irregular cuboctahedron 

In this section we investigate how the magnetic properties of the cubocta- 
hedron change if random variations of the exchange coupling parameters are 
introduced. This study is motivated by recent magnetization measurements 
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of the icosidodecahedral molecules {Mo72Fe3o} [H] and {Mo72Cr3o} [20] pub- 
lished in Ref. [28], which were interpreted by assuming random distributions 
of exchange parameters in a classical Heisenberg model description. 
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Fig. 6. Differential susceptibility as a function of applied field for the irregular 
cuboctahedron with s = 1/2. L.h.s.: dependence on the width A of the random 
distribution. Middle: dependence on the temperature kBT/\J\ = 0.25,0.5,0.75,1.0 
for A = 1.0|J[. R.h.s.: same as middle for A = 2.0|J|. 

We introduce variations of the exchange parameters of the Hamiltonian ([T]) 
by replacing the common nearest neighbor exchange parameter = J with 
values of a fiat random distribution J — 0.5A < J^^ < </ + 0.5A. Thus the 
mean exchange parameter is kept to be J. In order to gain sufficient statistical 
certainty we use ensembles of 10,000 spectra for realizations of the irregular 
cuboctahedron with s = 1/2; the results do not deviate from those for 1,000 
realizations. For larger s the production of sufficiently large ensembles is hin- 
dered by prohibitively many diagonalizations of larger matrices. 

Figure [6] shows the differential susceptibility that results from averages using 
distributions with various A. The figure on the l.h.s. compares dAi/dB at the 
rather low temperature of ksT = 0.1\J\ for A = 0, i.e. the regular cubocta- 
hedron, with A = 0.5|J|, A = 1.0|J|, and A = 2.0|J|. One clearly sees that 
the pattern which mainly originates from ground state level crossings does not 
change much for A = 0.5 |J| and A = 1.0 |J|. It needs a variation as large as 
A = 2.0| J|, i.e. ferromagnetic interactions occur, to qualitatively change the 
differential susceptibility function. The reason is that smaller variations do 
no alter the structure of low-lying energy gaps. The singlet-triplet gap, which 
is approximately 0.765|J|, does not vary very much for the ensembles with 
smaller A, and so does the singlet-triplet crossing which is determined by the 
singlet-triplet gap. It needs an appreciable variance of the exchange parameter 
distribution in order to impose large variations of the level crossing fields. 

The middle and the r.h.s. of Fig.[6]display dAi/dB for temperatures ksT /\J\ = 
0.25, 0.5, 0.75, 1.0 and A = 1.0| J| and A = 2.0| J|, respectively. As already ex- 
plained, there is only very little difference between the behavior of an irregular 
cuboctahedron with A = 1.0| J| (middle) and the regular one. For A = 2.0| J| 
(r.h.s.) the differential susceptibility is much more smeared out which includes 
an appreciable broadening at the saturation field. Considering the irregular 
cuboctahedron we can conclude that the magnetic properties are rather stable 
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against random fluctuations of tlie excliange parameters. Tliis means tfiat tfie 
striking beliavior especially of the experimental differential susceptibility of 
{Mo72Fe3o} and {Mo72Cr3o} which shows no signs of level crossings at all [2B] 
needs further theoretical exploration of the microscopic origin. 
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